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Abstract: We study the many-body theory of graphene Dirac quasiparticles interacting 
via the long-range Coulomb potential, taking as a starting point the ladder approximation 
to different vertex functions. We test in this way the low-energy behavior of the electron 
system beyond the simple logarithmic dependence of electronic correlators on the high- 
energy cutoff, which is characteristic of the large-A^ approximation. We show that the 
graphene many-body theory is perfectly renormalizable in the ladder approximation, as 
all higher powers in the cutoff dependence can be absorbed into the redefinition of a finite 
number of parameters (namely, the Fermi velocity and the weight of the fields) that remain 
free of infrared divergences even at the charge neutrality point. We illustrate this fact 
in the case of the vertex for the current density, where a complete cancellation between 
the cutoff dependences of vertex and electron self-energy corrections becomes crucial for 
the preservation of the gauge invariance of the theory. The other potentially divergent 
vertex corresponds to the staggered (sublattice odd) charge density, which is made cutoff 
independent by a redefinition in the scale of the density operator. This allows to compute 
a well-defined, scale invariant anomalous dimension to all orders in the ladder series, which 
becomes singular at a value of the interaction strength marking the onset of chiral symmetry 
breaking (and gap opening) in the Dirac field theory. The critical coupling we obtain in this 
way matches with great accuracy the value found with a quite different method, based on 
the resolution of the gap equation, thus reassuring the predictability of our renormalization 
approach. 
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1 Introduction 



The discovery of graphene, the two-dimensional material made of a one-atom-thick car- 
bon layer [1], has opened new possibilities to investigate fundamental physics as well for 
devising technological applications. The electron system has relativistic-like invariance 
at low-energies, mimicking the behavior of Dirac fermions in two spatial dimensions [2-4]. 
Moreover, the Coulomb repulsion between electrons constitutes the dominant interaction 
in the graphehe layer. This makes the low-energy theory to be a variant of Quantum 
Electrodynamics, but placed in the strong coupling regime as the ratio of to the Fermi 
velocity vp the electrons is nominally larger than one. 

There have been already several proposals to observe unconventional signatures of the 
interacting electrons in graphene. It has been for instance remarked that the interaction 
with impurities carrying a sufficiently large charge should result in anomalous screening 
properties of the graphene system[5-8]. Furthermore, it was also found long ago that 
the own e-e interaction in the layer should lead to a linear dependence on energy of the 
quasiparticle decay rate[9], as a consequence of the vanishing density of states at the charge 
neutrality point, and in agreement with measurements carried out in graphite[10]. 

More precisely, it has been shown that the graphene electron system has the proper- 
ties of a renormalizable quantum field theory, where the parameters flow with the energy 
scale[ll]. In this framework, it was implied that the Fermi velocity should scale logarith- 
mically towards larger values in the low-energy limit, what appears to be confirmed by 
recent experimental observations in graphene[12]. 

The graphene electron system is actually an example of electron liquid with strong 
many-body corrections, which depend significantly on the energy of the interaction pro- 
cesses. In practice, this is manifest in the logarithmic dependence on the high-energy 
cutoff needed to regularize the contributions to different quantities like the Fermi velocity 
or the weight of electron quasiparticles. In this type of electron liquid, one has to make 
sure that these divergences amount to the redefinition of a finite number of parameters in 
the system. In the context of quantum field theory, this property of renormalizability is 
crucial to guarantee the predictability of the theory as quantum corrections are taken into 
account. Otherwise, there is the possibility that the singular dependences on the cutoff 
cannot be absorbed into the redefinition of a finite number of local operators of the bare 
theory. This may happen when they take for instance the form of momentum-dependent 
log(|p|) corrections to local operators in the effective action, being then the refiection that 
the effective low-energy theory cannot be captured in terms of the local fields present in 
the original model. 

At this point, the best evidence of the renormalizability of the model of Dirac fermions 
in graphene comes from the study of the theory in the limit of large number N of fermion 
flavors [13], equivalent to the random-phase approximation (RPA). In this regime, it has 
been shown that all the cutoff dependences of the theory can be absorbed into redefinitions 
of the Fermi velocity and the weight of the electron quasiparticles, to all orders of the 
perturbative expansion in e'^/vp [13] (for other studies of the 1/A^ expansion in graphene, 
see also Refs. [14] and [15]). Anyhow, many-body corrections only exhibit at large N a 
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simple logarithmic dependence on the energy cutoff, which makes rather straightforward 
the renormalization of the theory at this stage. 

In this paper we adopt an approach that is opposite in many aspects to that of the 
large- expansion, and that is able to probe the structure of the many-body corrections 
with arbitrary large powers of the cutoff dependence. That is based on the sum of the 
series of ladder diagrams, that we apply to different interaction vertices of the theory. 
Within this approach, we will be able to show that the divergent dependences on the cutoff 
can be reabsorbed in a finite number of parameters of the interacting theory, including 
the renormalization of the scale of different bilinears of the Dirac fermion fields. These 
renormalized quantities will prove to be independent of any infrared scale (Fermi energy, 
external momenta), making the low-energy limit of the many-body theory perfectly well- 
defined even at the charge neutrality point. 

From a practical point of view, the motivation for focusing on the sum of the ladder 
series lies in that it encodes the most divergent diagrams at each level of the perturbative 
expansion for the undoped electron system[16]. This makes highly nontrivial the process of 
renormalization, by which one has to remove in general divergent corrections that behave 
like the n-th power of the logarithm of the cutoff, when looking at the n-th perturbative 
level. In practice, we will illustrate the usefulness of the renormalization approach in 
the computation of observables like the anomalous dimensions of composite operators, 
which become determined just by the value of the renormalized coupling constant. This 
will allow us to address in particular the question of the dynamical breakdown of the 
chiral symmetry in the electron system [17-29], which can be characterized in terms of the 
singular behavior of the corresponding anomalous dimension at a certain critical value of 
the coupling constant[30, 31]. 

2 Dirac many-body theory 

Graphene is a 2D crystal of carbon atoms forming a honeycomb lattice, such that its low- 
energy electron quasiparticles are disposed into conical conduction and valence bands that 
touch at the six corners of the Brillouin zone[4]. Of all six Fermi points, there are only two 
independent classes of electronic excitations. Thus, the low-energy electronic states can 
be encoded into a set of four-dimensional Dirac spinors {il^i}, which are characterized by 
having linear energy-momentum dispersion e(p) = Wi?|p|. The index i accounts for the two 
spin degrees of freedom, but may also allow to extend formally the analysis for a higher 
number N of Dirac spinors. The kinetic term of the hamiltonian in this low-energy theory 
is given by 



where t/^^ = ipj'yo and {70-} is a collection of four-dimensional matrices such that {7^, 7;/} = 
2 diag(l, — 1, — 1). They can be conveniently represented in terms of Pauli matrices as 
70,1,2 = (<73; CTsCTi, C730"2) ® c^s, where the first factor acts on the two sublattice components 
of the honeycomb lattice and the second factor operates on the set of two independent 
Fermi points. 




(2.1) 
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In this paper we focus on the effects of the long-range Coulomb interaction in the 
graphene electron system. The density of states vanishes at the Fermi points connecting 
the conduction and valence bands, so that a sensible starting point for the e-e interaction 
is given by the unscreened potential V{r) = e^/47r|r|. The long-range Coulomb repulsion 
governs actually the properties of the electron system at low energies, since it is the only 
interaction that is not suppressed, at the classical level, when scaling the many-body theory 
in the limit of very large distances. If we add to (2.1) the contribution from the Coulomb 
interaction, we get the expression of the full hamiltonian 

H = -ivF [ d\ V^i(r)7 • VV'i(r) + ^ / d\i [ dVs p{ri)- ^ -p(r2) (2.2) 

J 8-71 J J |ri-r2| 

with p{r) = ipj^{r)^o'tpi{r). The total action of the system is 

S = j dt j d^r^i{r){i-iQdt + ivFl •V)^i{T) 

-^JdtJ d^n I d\2 p(ri)^^;-^/5(r2) (2.3) 

This action is invariant under the combined transformation of the space and time variables 
and the scale of the fields 

t' = st , r' = sr , = (2.4) 

This means in particular that the strength of the interaction is not diminished (contrary 
to the case of a short-range interaction) when zooming into the low-energy limit s — t- oo. 

This analysis shows that the Coulomb repulsion mediated by the long-distance l/|r| 
potential is the only interaction that may prevail in the low-energy regime of the electron 
system. It is clear that any other e-e interaction without the l/|r| tail, as those that arise 
effectively from phonon exchange, will be suppressed at least by a power of 1/s under 
the change of variables (2.4). The Dirac field theory with long-range Coulomb interaction 
has indeed the property of being scale invariant at this classical level, which provides a 
good starting point to investigate the behavior of the many-body corrections upon scale 
transformations towards the long- wavelength limit s — )■ oo. 

In fact, the many-body theory does not preserve in general the scale invariance of the 
classical action (2.3), as a high-energy cutoff A has to be introduced to obtain finite results 
in the computation of many-body corrections to different observables. The analysis of the 
cutoff dependence of the many-body theory provides deeper insight into the effective low- 
energy theory. If the theory is renormalizable, it must be possible to absorb all powers of 
the cutoff dependence into a redefinition of the parameters in the action (2.3). This should 
be therefore modified to read 

S = Zkin j dtj d^r i^i{r){i-fodt + iZ^ vpl ■ V)V'i(r) 
-Zint|^ j dt j d\i j d\2 p{ri)j-^—^p{r2) (2.5) 
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Figure 1. Diagrammatic representation of the ladder approximation for tlie electron self-energy. 



The assumption is that Z^in, and Zint (and other renormalization factors for composite 
operators that do not appear in (2.3)) can only depend on the cutoff, while they must be 
precisely chosen to render all electronic correlators cutoff independent. 

The renormalizability of the graphene Dirac field theory is a nontrivial statement, 
since it amounts to the fact that that all the many-body corrections depending on the 
high-energy cutoff must reproduce the structure of the simple local operators that appear 
in (2.3). It turns out for instance that many of the individual contributions to a given 
correlator have dependences in momentum space of the form log"(|p|) log™' (A). These are 
nonlocal corrections that cannot be reabsorbed into the action (2.5), and the fact that all 
these nonlocal terms cancel out in the final result for a correlator is a remarkable property 
of a renormalizable theory. Non-renormalizable theories have in this regard an essential 
lack of predictability, as log"(|p|) corrections do not make viable the characterization of 
the low-energy effective theory in terms of a few local operators, which may be in turn the 
reflection that it is not actually captured by the original fields formulated in the many-body 
theory. 



3 Electron self-energy and Fermi velocity renormalization 

We first consider the cutoff dependence of the electron self-energy in the ladder approxima- 
tion. We define this approach in terms of the self-consistent equation represented in Fig. 
1. Diagrammatically, it corresponds to build the electron self-energy by iteration in the 
number of "rainbow-like" interactions between the Dirac fermion lines. A similar approach 
will be used afterwards to define the ladder approximation for the vertices of the charge 
and current density operators. 

Before dealing with the actual ladder series, we establish our representation of the free 
propagators by describing the computation of the lowest-order self-energy diagram. The 
free propagation of the Dirac fermions corresponds for instance to the expectation value 

(V'i(k,a;)'0j(k,a;))fi.cc = iGo(k,a;) 

. -7ow + VFl • k 
= I 5 n-p5 — (3.1) 

—ixl'^ + VpVL'^ — IT] 

On the other hand, the interaction lines stand in momentum space for the product of —i 
times the Fourier transform of the Coulomb potential, that turns out to be in two spatial 
dimensions 

v^(q) = ^ (3.2) 
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The first-order electron self-energy diagram, that we will denote by Si(k), needs to 
be regularized by introducing a high-energy cutoff A in the momentum integrals. We have 
actually 

2 J [ZTry J ZTT —UJ^ + VpP"^ — IT] |k — p| 

which leads to a contribution proportional to 7 • k that must be logarithmically divergent 
by simple dimensional counting. If we bound the integration in momentum space such that 
vf\p\ < A, we get the result 

Si(k)«-^7-k logA (3.4) 

that corresponds to the well-known renormalization of the Fermi velocity by the Coulomb 
interaction in the Dirac many-body theory[ll, 13]. 

From now on, we will choose a convenient regularization method to compute the di- 
vergent as well as the finite corrections to electronic correlators at each perturbative level. 
That consists in the analytic continuation in the number of space dimensions [32], by which 
the momentum integrals are computed at dimension D = 2 — e [11]. With this method, 
dependences on log A are traded in general by 1/e poles. In the above instance of the 
electron self-energy, we get after integration over ujp 

zS,(k)=4 / T^^T-Pn^ (3-5) 
4 J (2tt)^ IpI |k- p| 

where eo is a parameter whose dimensions are given by an auxiliary momentum scale ^ 
through the relation 

eo = /i^/2e (3.6) 
The calculation then proceeds as follows: 

Si(k) = -jz / TT^f-P / dx- 



AttJ {2tt)D ' "^Jo (k-p)2x + p2(l -x) 

d^p , f\ xi/2(l_:c)-V2 
7 • k / dx- 



An J (27r)^ ' Jq p2 + k'^x{l - x 
4.. u f'^ Til-D/2) 



1^7 • k / dx 



(47r)^/2 (k2x(l - x))i-D/2 

_ ^ (4vr)-/2r(ie)r(3^)r(i^) 

" (47r)2^"^ |k|^ r(2-6) ^"^-^^ 

From the latter expression we find the pole as e — )• 

^.W»I^7-ki (3,8) 

We are anyhow interested in the result of computing the electron self-energy in the 
ladder approximation defined in Fig. 1. It is easily realized that the solution Siadder(k) of 
the self-consistent equation must have the structure 

Sladdcr(k) =/(k)7-k (3.9) 
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with a scalar function /(k). The self-energy then satisfies 

i5]laddcr(k) = ii;i(k) + ^ J I ^Sladdcr(p) ^ , .9 ^ ^'.2 ^"9" ,■„^9 n. ^ „| (3-10) 



The solution of Eq. (3.10) reflects a particular feature of the graphene many-body 
theory in the static limit (i.e. when the effective interaction is supposed to be frequency 
independent). It can be checked that the second term in the self-consistent equation iden- 
tically vanishes, as a result of performing the integration over the frequency variable, and 
irrespective of the actual momentum dependence of Siaddcr(k). We have indeed, by per- 
forming a Wick rotation to imaginary frequency uJk = ~i^k^ 



I 



/ ^\_2 2 9X9 =Q 3.11 



This result implies that the solution of Eq. (3.10) must coincide with the first-order con- 
tribution 

Sladdcr(k) = Si(k) (3.12) 

This vanishing of higher-order corrections to the electron self-energy in the ladder 
approximation can be actually seen as the consequence of a wider symmetry operating in 
the graphene many-body theory. We can extend the sum of self-energy diagrams in the 
ladder series to include contributions where the electron lines in the ladder diagrams are 
corrected by the own electron self-energy. This leads to the sum of a much broader class of 
diagrams, that are encoded in the self-consistent exchange approximation depicted in Fig. 
2. If we represent the electron self-energy in this approach by 

SscEx(k) = /(k)7-k (3.13) 

the self-consistent equation can be written as 

zLscEx(k) = -— / — / 17-70 , .,3 g -TOTT — -| (3.14) 

2 J [2tT)^ J 2tT -oj^ + (yp + /(p))2p2 - IT] |k - p| 

The key observation is that, for the same reason that the final expression for the first- 
order self-energy (3.3) does not depend on the Fermi velocity vp, the integral in Eq. (3.14) 
turns out to be independent of the function /(k). One can for instance redefine the scale 
of the frequency from ujp to (1 + f{p)/vF)ojp, in such a way that the integrand in Eq. 
(3.14) falls into the corresponding first-order expression in (3.3). This proves that, also in 
the more comprehensive self-consistent exchange approximation, the electron self-energy 
coincides with the first-order result 

SscEx(k) = Si(k) (3.15) 

In its simplicity, the result expressed in Eq. (3.15) accounts for the vanishing of a 
vast class of corrections to the electron self-energy in graphene. It can be interpreted as a 
kind of no-renormalization theorem that protects the Fermi velocity from being modified 
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Figure 2. Diagrammatic representation of the self-consistent exchange approximation for the 
electron self-energy. 
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Figure 3. Diagrammatic representation of the vertices for (a) total charge density, (b) current 
density, and (c) staggered charge density. 

by higher-order effects, which remains valid under the assumption of static (frequency- 
independent) screening of the Coulomb interaction. As we will see, this translates into 
a remarkable cancellation of corrections to the vertex of the current density operator, in 
a nontrivial check of a gauge invariance that is hidden in the original formulation of the 
theory. 

4 Charge and current density correlations 

We study next the way in which many-body corrections dress the different interaction 
vertices in graphene. This includes the inspection of the own Coulomb interaction, that we 
analyze by looking for corrections of the coupling to the total charge density. The Dirac 
field theory allows anyhow for the consideration of more general vertices that take into 
account the pseudospin current 7 and the spinor structure of the fermion fields. We will 
pay attention in what follows to the vertices for the total charge, the pseudospin current, 
and the staggered (sublattice odd) charge density, which are represented in Fig. 3. 

The analysis of the vertices for the total charge and the pseudospin current is particu- 
larly relevant since, if we think of them as operators that may be switched on in the action 
of the electron system, it becomes clear that they should be related by gauge invariance to 
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Figure 4. Self-consistent diagrammatic equation for a generic vertex Ti in the ladder approxima- 
tion. 

the terms in the kinetic action. That is, we can start with an extended action given by 



where and A play the role of auxiliary fields mediating the interactions of the total 
charge and the pseudospin current. A gauge transformation of the Dirac fields 



amounts to a shift of the auxiliary fields Aq and A. Thus, the invariance of the many-body 
theory under (4.2) can be tested by checking that the renormalization factors Z^^^ and 
match with the respective factors from the renormalization of the electron self-energy. 

This question of the gauge invariance is an interesting point regarding the graphene 
many-body theory, as it was shown long ago that the four-fermion interactions in the 
graphene electron system can be obtained from a suitable projection of the full relativistic 
interaction mediated by photons in three spatial dimensions [11]. It has been actually proven 
that the renormalization of the theory, when carried out to first order in perturbation 
theory, is consistent with the above mentioned gauge invariance. In the present instance, 
we will also use the renormalization properties of the vertices to check the underlying 
gauge invariance to higher orders in the ladder approximation supplemented by electron 
self-energy corrections. 

As in the case of the electron self-energy, we define the ladder approach for the vertices 
by means of a self-consistent equation, represented now in Fig. 4. In principle, one can 
solve the equation by means of an iterative procedure, ending up with the equivalent of a 
ladder series for the different vertices. We will then improve this diagrammatic approach in 
a second stage, by assuming that the internal fermion lines in the self-consistent equation 
are themselves corrected by the electron self-energy, which will prove to be crucial to 
preserve the gauge invariance of the theory. 

4.1 Charge density vertex 

We define the vertex for the total charge density in frequency and momentum space as 






(4.2) 




(4.3) 
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where p is given in real space by 



p(r) = V'i(r)7oV'i(r) (4.4) 

and the right-hand-side of (4.3) is computed by considering only the one-particle-irreducible 
vertex diagrams. In this way, the possible renormalization required to render Vq cutoff 
independent should amount to a simple multiplication of the vertex by a factor, that is the 
same Z[^^ appearing in (4.1). 

The vertex Vq is a dimensionless quantity, which means that, in order to isolate the 
singular dependence on the cutoff, it is enough to study the limit q — t- and ujq — t- 0. Then, 
the self-consistent equation represented in Fig. 4 becomes 

r„{0,0;.,.0 =70 + 4 / ^|?r,,O,0;p,.,)^-;^ii^|^^ ,4.5, 

It is clear that the solution of (4.5) cannot depend on the frequency ook of the external 
fermion lines. Therefore, the integral at the right-hand-side of the equation must be iden- 
tically zero, for the same reason that the integral in Eq. (3.10) was also vanishing. This 
means that the vertex Fq is independent of the cutoff in the ladder approximation. Repeat- 
ing here the argument at the end of the last section, it turns out that the same statement 
holds true even when the fermion propagators in (4.5) are corrected with the electron self- 
energy (3.7). Again, the integral in (4.5) vanishes irrespective of the momentum-dependent 
corrections to vp, leaving Fq cutoff independent in this approach. 

The cutoff independence of Fq agrees with the absence of wavefunction renormalization 
(Zkin = 1) in the self-consistent exchange approximation applied to the electron self-energy. 
The trivial result 

Zkin = ^(„t = 1 (4.6) 

is the first check of the gauge invariance of the theory. The vanishing of the many-body 
corrections lies in this instance in the particular structure of the ladder approximation and, 
in this regard, it is a result that holds even after dressing the interaction with the static 
(frequency-independent) RPA screening of the Coulomb potential. 

4.2 Current density vertex 

The irreducible vertex for the current density is defined in this case by 

rc(q, ojq\ k, Wfc) = (Pc(qi Wg)?/'i(k + q, Wfc + a;g)V'j(k, a;fc))ipi (4.7) 

where the current density operator is given in real space by 

p,(r) =V^,(r)7^i(r) (4.8) 

We anticipate the fact that the computation of the vertex may give rise to dependences 
on the high-energy cutoff, that are supposed to be absorbed in the renormalization factor 
7" 
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The vertex Tc is a two-dimensional vector, but its analysis can be greatly simplified by 
considering as before the limit q — t- and Uq — t- 0. The self-consistent equation depicted 
in Fig. 4 takes then the form 

re(O,0;k,Wfc) = 

We resort at this point to an iterative resolution of (4.9), by which we can obtain 
a recursion between consecutive orders in the expansion of the vertex in powers of the 
interaction strength. This procedure shows that Tc{O,0;'k,u)k) has a part proportional to 
7 and another contribution proportional to k(7-k). Prom dimensional arguments, one can 
see that the solution of (4.9) must take the form 

(OO \ OO / 

1 + E ^0 7^7 + • Hk) Yl ^0 7^^ (4-10) 
n=l ' ' / n=l ' ' 

where we have called nk = k/|k| and 

If we insert for instance a given order of the expansion with coefficient r„ inside the 
integral in Eq. (4.9), we get 

-ijoZJp + VFJ ■ p r„ -ijoZJp + VF-f ■ p 1 



el f d^, 
2 J (27r; 

el 


p dLjp 
)D 2. ^° 

{ d^'p 


el 


f d^p 




i2Tr)D ^ 



(^^ + w^p2 '''IpI"^ i^^ + ?;|,p2 ^°|k-p| 



p|3+ney |k _ p| 
( _ ( ^^ r(l + f) ^-1/2(1 _^)-(l-n.)/2 

(7 "piT-npjj^^^i^^y^ ((k-p)2x + p2(l-x))i+W2 

(4vr)^/2 /r(i^e)r(i^)r(i^) r (M:i.) r (1^) r (^) 

7 -^0 



k|("+i)^ 4 I (i±^^) r (i - 2y?r (3±i^) r (2 



(n+2)e 
2 



r (1 + isfij) r (tMili) r , , 



On the other hand, by inserting any term of the expansion (4.10) with coefficient inside 
the integral of the self-consistent equation, we get always a vanishing result due to Eq. 
(3.11). We obtain therefore the recurrence relations 



rn+l 



(4vr)^/2 /r (H±ie) r (i^) r C-^) r m^e) r (i^) r (3^) ' 



V^T {^) r (i - 2V^r {^) r (2 



(n+2)£ 
2 



(4.13) 
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(4.r/2r(i + ^e)r 5^ r(i-) 



v^r (3^) r (2 



(n+2)e 
2 



We observe from (4.13) that the expansion of the vertex develops increasing divergences 
in the cutoff, that manifest as poles in the limit e — t- 0. The key point is whether these 
divergences can be absorbed by a suitable renormalization factor. The current density Pc(r) 
is not an elementary field of the many-body theory, which means that its correlators need 
to be renormalized by appropriate rescaling of the own current density. Alternatively, if we 
include the composite field with its own coupling in the action (4.1), it is the renormalization 
factor Z[^^ which needs to be adjusted to render the correlators cutoff independent. Then, a 
renormalized vertex Tc^rem finite in the limit e — )• 0, has to be obtained by the multiplicative 
renormalization 

r,,ren = Z!'^,T, (4.15) 
The renormalization factor may have in general the structure 

^1 = 1 + 1;^ (4-16) 

1=1 

in terms of the dimensionless physical coupling 

Attvf 



A = (4.17) 



It is a nontrivial fact that all the poles in Tc may be canceled against multiplication by 
Z[^^, allowing only for the dependence of the coefficients Cj on the coupling constant. We 
have checked that this is indeed the case, up to the order A^^ we have been able to carry 
out the numerical computation of the vertex. We have found for instance for the first terms 
in the expansion (4.16) 

ci(A) = -^A - ^(1 + log(16))A2 - ^ (-1 + 31og2(4) + log(64)) A^ 

-9 - 721og(2) + 3841og2(2) + 1281og2(2) log(16) + 12C(3) 4 

A 



24576 

3 - 301og(2) - 601og^(2) + 4001og^(2) + 4001og'^(2) + 6C(3) + 61og(16)C(3) 

24576 

+ ... (4.18) 

-X^ + ... (4.19) 



13 + 641og(2) - 5281og2(2) - 176 log2(2) log(16) - 12C(3) 5 



98304 

1 3 l + 41og(2) 4 -5 + 721og(2) + 1441og^(2) 5 

C3(A) = -—A - 2048 ^ 98304 ^ + " " " ^^'^^^ 

C4(A) = ^A^ + i±lM^A^ + ... (4.21) 

^ ^ 6144 24576 ^ ^ 

C5(A) = — A^ + ... (4.22) 

^ 122880 ^ ^ 
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The important point about this result for Zj"^ is that it does not depend on the momenta 
of the vertex Tc- This means that it represents a local divergence as e — ?• 0, and it can be 
therefore understood as the renormalization of a local operator in the action (4.1). 

Another important consequence of the actual expression of the functions Cj(A) is that 
observable quantities derived from Z^^^, as for instance the anomalous dimension 7c of the 
current operator, turn out to be finite in the limit e — )• 0. The dimension 7c measures in 
particular the anomalous scaling of under changes in the units of energy and momentum 
in the system. The vertex Tc is formally a dimensionless quantity, but the renormalization 
process introduces scale dependence on the auxiliary momentum /x, in such a way that 

Tcren ~ (4.23) 

The anomalous scaling of the vertex comes only from the dependence of the renormalization 
factor Zj'^^ on the fi scale, so that 

Assuming the general structure (4.16), it is in general a nontrivial fact that the anoma- 
lous dimension 7c computed from Z'(^^ may become finite in the limit e — )• 0. The depen- 
dence on the scale /i is encoded in the equation 

= -A (4.25) 

We can then express Eq. (4.24) in the form 

Zl,di, dX 

mt j=o 



7c 



Alternatively, we can write the above equation as 



c,; \ , \ - dCi+i 

i=l ^ / i=0 



Assuming the finiteness of the anomalous dimension in the limit e — )• 0, we get [32] 

.. = -A^| (4.28) 

and the consistency conditions for the cancellation of the poles at e = 

If -It 

Quite remarkably, it can be seen that the expressions in (4.18-4.22) satisfy identically 
the conditions (4.29). This holds for the functions Cj(A) which we have obtained analyt- 
ically up to order A^. We have been also able to compute numerically their power series 
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expansion up to order A^®, checking that the equations (4.29) are verified order by order 
with the precision allowed by the calculation. This provides a very strong evidence of the 
renormalizability of the theory, implying that observable quantities like 7c can be computed 
from renormalized correlators to obtain finite results, dependent only on the value of the 
physical coupling constant. 

On the other hand, we note that the result for Fcrcn is drastically modified when the 
electron self-energy corrections are included in the calculation of the vertex. As we have 
already seen in Sec. 3, the main effect of the electron self-energy is to renormalize the value 
of the Fermi velocity vp. At the level discussed in that section, the self-energy corrections 
amount to perform the replacement in the inverse of the Dirac propagator 

jQUJ — vp-f-p — )• 'Jquj — vpj • p — T,i{p) (4.30) 

with Si(p) given by Eq. (3.7). It is then clear that the electron self-energy diagrams can 
be incorporated to the ladder approximation encoded in Eq. (4.9) simply by trading the 
constant vp hy an effective Fermi velocity 

MP) = v. + S^^-r^'^-^^^^w^W ^'-''^ 

The replacement of Ui? by vp{p) in the formulas has the effect of iterating in the number 
of self-energy diagrams inserted in the electron and hole propagators building the vertex. 
The electron self-energy corrections contribute therefore to supplement the ladder series 
previously considered, greatly improving the diagrammatic approach for the vertex Tc- 

In order to compute the renormalization factor Zj'^^. in the ladder approximation with 
effective Fermi velocity vp{p), it is convenient to expand the factor l/vp(p) in powers 
of Cq inside the integral of Eq. (4.12). The vertex Tc still admits a solution like that in 
Eq. (4.10), as each order in the perturbative expansion can be represented in terms of the 
precedent by integrals of the type shown in Eq. (4.12). The power series in Cg contains 
now more poles in the e parameter, as a result of the divergent behavior of vp{p) in Eq. 
(4.31). The poles coming from the electron self-energy corrections can be however removed 
at once by the renormalization of the Fermi velocity 

Vp = ZyVp^ren (4.32) 

with 

Z^ = l + bi- (4.33) 
e 

The coefficient needs simply to be adjusted to 

61 = -—^ (4.34) 

lOVrUF.ren 

leading then to a finite vf{p) written in terms of vp^^cn- 

We have again a general structure for the renormalization factor in this improved 
approach 

ZL = l + f:~^ (4.35) 

i=l 
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The remarkable result is that, after writing the perturbative expansion for the vertex as a 
power series in the renormalized coupling 

A = (4.36) 

one needs just a simple first-order 1/e term in (4.35) to get rid of all the poles in V^- That 
is, the renormalized vertex becomes finite in the limit e — )• with the choice 

ci(A) = -^A (4.37) 
q(A) = i>2 (4.38) 

The simple pole structure of Zj"^ in the improved ladder approximation implies the 
result 

= Zl, (4.39) 

We may see in this relation a nontrivial link between the renormalization of the ivp'f ■ V 
kinetic term and that of 7 • A in the action (4.1). This feature points at a symmetry 
that is characteristic of a gauge invariant theory, and it can be explained in our case from 
inspection of the diagrams contributing to the self-energy discussed in Sec. 3 and to the 
vertices computed in this section, as we show next. 

4.3 No-renormalization of charge and current density operators 

The result (4.39) implies, together with (4.6), the preservation of the gauge invariance in 
the renormalized action (4.1). In this respect, there are actually Ward identities that can 
be derived from general principles, relating different vertex functions of the theory. This 
is stressed for instance in Ref. [33], where it has been also emphasized the suitability of 
the dimensional regularization method to preserve the gauge symmetry of the theory. We 
show here that the two relevant identities between the electron self-energy and the vertices 
Fq and Tc can be easily obtained in the framework of the present many-body approach. 

The main idea is that any electron self-energy correction can be converted into a 
contribution to the vertex Tq or Tc by taking the derivative with respect to the external 
frequency uj^ or the external momentum k of the self-energy. This fact relies on the 
expression of the derivatives of the free Dirac propagator 

d 1 



duJk 70 (^fc - Wp) - vp^ ■ (k 



^ 70 r^ V (4-40) 



7o(wfc - ujp) - VFl • (k - p) 7o(wfc - Wp) - VFl • (k - p 
1 d 1 



VF 9k 7o(a;fc - ujp) - vfI • (k - p) 

= ^- y i (4.41) 

7o(wfc - Wp) - VFl • (k - p) 7o(wfc - ^p) - VFl • (k - p) 

In any correction to the self-energy S(k, w^), one can choose the external frequency and 
momentum to circulate along the fermion lines that connect the two outer vertices of the 
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Figure 5. Schematic representation of tiie Ward identity between the electron self-energy and the 
current density vertex in the ladder approximation. 

diagram. Thus, taking the derivative with respect to or k implies cutting any of those 
internal lines in two pieces and inserting a vertex with the respective 70 or 7 matrices. 

The above construction becomes clear if one has in mind the diagrams building the self- 
consistent exchange approximation considered in Sec. 3. Taking the derivative of SscEx(k) 
with respect to uj^ gives identically zero, which is consistent with the absence of corrections 
to Tq in the ladder approximation (at zero momentum transfer). On the other hand, the 
derivative with respect to k gives rise to two different types of diagrams contributing to 
Fc, as shown in Fig. 5. Part of them corresponds to the kind of contributions that we were 
considering in the ladder approximation to Fc without electron self-energy corrections, as 
seen in the upper right of the figure. But the other part of the diagrams consists of vertex 
corrections with self-energy insertions in the electron and hole internal lines, as illustrated 
in the lower right of Fig. 5. This shows that the differentiation of SgcEx(k) generates 
actually the whole set of vertex corrections in the ladder approximation supplemented 
with electron self-energy corrections. 

We can then write a Ward identity of the form 



This identity implies that the renormalization of the vertex Fc is dictated by that of the 
Fermi velocity vp- In this regard, the result = Zj'^^ found above becomes a natural 
consequence of Eq. (4.42). Alternatively, these findings also stress the fact that the electron 
self-energy corrections cannot be neglected in a consistent approximation to the many-body 
theory of graphene, as they play a crucial role to build a gauge invariant effective action 
with the structure given by Eq. (4.1). 




(4.42) 



5 Staggered (sublattice odd) charge density correlations 



m 



We may also consider the renormalization of the staggered charge density operator anti- 
symmetric under the exchange of the two sublattices of the graphene honeycomb lattice 

P3(r)=?,(r)Vi(r) (5.1) 

We will define the corresponding vertex by 

r3(q,tjg;k,a;fc) = {p3{ci,ujg)ipi(k + <i,ujk + ujg)Tpi(k,ujk))iPl (5.2) 

where IPI denotes again that we take the irreducible part of the correlator. 

The vertex has a clear physical significance as it enters in the correlations of the 
staggered charge psir), which is the order parameter for chiral symmetry breaking in the 
many-body theory. A nonvanishing expectation value (p3(r)) 7^ is the signal that a mass 
is dynamically generated for the Dirac fermions. This means that their hamiltonian gets 
effectively a term of the form 

: I (fri^i{r)^IJi{r) (5.3) 

With the mass term, the conduction and valence bands loose the perfect conical shape 
about the charge neutrality point, and a gap opens in the electronic spectrum. This trend 
of symmetry breaking is similar to that discussed long ago in Quantum Electrodynamics 
in two spatial dimensions [34-38]. In the present context, the dynamical mass generation is 
also driven by the interaction, in such a way that the condensation of P3(r) may proceed 
depending on the value of the coupling /A.ttvf (and also on the number of fermion flavors, 
in a theory with a number N of different fermion species). 

5.1 Staggered charge density vertex in ladder approximation 

We deal first with the vertex r3 in the ladder approximation, which is given again by 
the self-consistent equation represented diagrammatically in Fig. 4. As we are mainly 
interested in the cutoff dependence of the vertex, we can take in particular momentum 
transfer q = and ujq = 0. Given that r3 must be anyhow proportional to the identity 
matrix, we have 

--jpUp + VFJ • P ^» n -joUJp + VF2^ r3(0,0;p,^p) 

^0 2,22 r3(0,0;p,Wp) 2,22 -70 = 2 I 2 2 ~ (^•'^) 

—iOp + VpP"^ — IT] —UJp + VpP^ — 17] —iOp + U^p — 11] 

Thus, the self-consistent equation for the vertex in the ladder approximation becomes 

p2 f fiD ] 1 1 
r3(0,0;k,^,,) = l-^^ / y^^r3(0,0;p,a;,) ^ ^ ^ . , (5.5) 

Eq. (5.5) can be further simplified by noticing that the solution cannot depend on the 
frequency oj^- We end up then with the equation 

r3(0,0;k,wfc) = l + ^ / -^r3(0,0;p,a;fc) ] . ..^ , (5.6) 
4 J {2tt)^ VF\p\\k — p\ 



-17- 



Prom dimensional arguments, the solution of Eq. (5.6) can be expressed in the form 



r3(0,0;k,a;,) = l + E^o^ (5-7) 



n=l ' ' 

with Ao = Cq/A-kvp- Each term in the series (5.7) can be obtained from the previous one, 
noticing that if we insert the general term in the integral at the right-hand-side of Eq. (5.6) 
we get 

" d^p 1 I 1 (4vr)-/^r(rr^e)r(i^)r(i^) 1 

(2^)^ IpI'^^vfIpI |k-p| ° 4 (i^) r (1 - ^e) |k|(™+i)^ 

(5.8) 

We observe that the result of the integral diverges in the limit e — t- 0, leading to a sequence 
of higher-order poles in the e parameter as we look at higher perturbative levels in the 
solution (5.7). 

The poles that appear in the computation of at e = must be reabsorbed by a 
suitable redefinition in the scale of the operator P3(r). This is a composite field, suscep- 
tible of being renormalized by a factor Zm which is independent of the renormalization 
of the elementary fields in the action [39]. This redefinition of /03(r) translates into the 
multiplicative renormalization of the vertex 

r3,rcn = ■^mr3 (5.9) 

The general structure of Zm must be 

^^ = i+f;^ (5.10) 

in order to absorb all the poles generated by the recurrence relation (5.8). 

A nontrivial check of the renormalizability of the theory is that the vertex Fs^ren must 
have a finite limit as e — t- 0, after making an appropriate choice of functions di{X) depending 
only on the coupling constant A. We have seen that this is the case up to the order A^^ 
we have pursued the numerical calculation of r3, finding a set of di{\) that do not depend 
on the momentum k of the vertex. The perturbative expansion of the functions can be 
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computed analytically with some effort up to order A®, leading to 

501og^(2) + 31og(2)C(3) ^5 _ 46081og5(2) + 4801og^(2)C(3) + 5C(5) 
96 5120 

.,,A) . i + i l„g(2) A» + A A* + 

^ ia21og''(2)+91og(2X(3) ^ 

*(^) - J (^) ^* - ^ '-^<^) - """11;^^'^' - ■ ■ (-3) 

di{\) = —x^ + — iog(2) + — log2(2) A^ + . . . (5.14) 
^ ^ 384 192 ^ 768 ^ ^ ^ ' 

An internal consistency check of the renormalizable theory is that, as in the case of the 
vertex Fc^reiD the computation of physical observables like the anomalous dimension 7^ of 
the operator has to provide a finite result in the limit e — t- 0. 7^ is defined in terms of 
the anomalous scaling of the vertex as a function of the dimensionful parameter ^ in the 
renormalized theory, 

r3,ren ~ /U^" (5.17) 

The dependence on ^ arises from the renormalization factor Z^, so that 

7™ = f (5.18) 
We can follow the same derivation as in Eqs. (4.25)-(4.27), ending up in the equation 

\ i=l / i=0 

From Eq. (5.19) we obtain a finite answer at e = for the anomalous dimension[32] 

7m = -A^di(A) (5.20) 
provided that the recurrence relations 

j^di+i{X) = di{X) ^di{\) (5.21) 

are identically satisfied. 

We have verified that the conditions (5.21) are fulfilled, up to the order A^^ we have 
computed numerically the perturbative expansion of the functions di{\). This means that 
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Figure 6. Plot of the absolute value of the coefficients dj"-* in the expansion of di(A) as a power 
series of the renormalized coupling A. 



the anomalous dimension 7^, is perfectly well-defined by Eq. (5.20) in the present frame- 
work. Regarding the actual calculation, we have found that the perturbative expansion 



di(A) 



n=l 



(5.22) 



behaves as a power series with a finite radius of convergence Ac- The geometric growth of 

(n) 

the coefficients d\ is illustrated in the plot of Fig. 6. The radius of convergence can be 



obtained by computing the ratio between consecutive orders of d\ 
converges towards a limit value d. It then turns out that 



(n) 



di{\) = ^(i"A" + regular terms 



and noticing that it 



(5.23) 



n=l 



An excellent fit of the n-dependence of d^-^^^^ /d^^^ is achieved by assuming the scaling 
behavior 



d- 



(n) 



d' d" d" 



(n-1) 



d+- + 



+ 



+ 



d 



n 



(5.24) 



(5.25) 



We obtain in this way an estimate of the radius of convergence 

Ac = i 0.456947 

d 

The singular behavior of the anomalous dimension 7^ is the manifestation of the 
divergence of the vertex Fa^ren at the critical coupling Ac. It also implies the divergence of 
the correlators of the staggered charge density P3{r) as these need to be multiplicatively 
renormalized by factors of Zm, which means that their anomalous dimensions are given 
by multiples of 7^. What we find therefore at the critical point Ac is the signature of the 
dynamical breakdown of the chiral symmetry in the Dirac theory. The predicted critical 
coupling in Eq. (5.25) turns out to be a very accurate approximation to the critical value 

87r2 



A. 



(5.26) 
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which has been obtained by Gamayun et al. in Ref. [27] by a quite different approach, 
consisting in the self-consistent resolution of the gap equation for the Dirac fermions (for 
another connection between the branching point of the gap equation and the singularity 
in the vertex, see also Ref. [40]). This remarkable coincidence between the results of two 
completely different methodologies can be taken as the reflection that they are encoding 
at the end an equivalent sum of many-body corrections, thus providing a nice check of the 
reliability of our computational framework based on the ladder approximation to vertex 
functions. 

5.2 Staggered charge density vertex supplemented with self-energy correc- 
tions 

We have anyhow to keep in mind that the electron self-energy corrections need to be 
incorporated to reach sensible results for rs^i-eru as we learned from the renormalization 
of the current density vertex. In the present case, the ladder approximation can be also 
improved by inserting the series of electron self-energy diagrams in the internal electron 
and hole states of the vertex. Prom a computational point of view, this can be achieved 
by replacing the constant vp in the integrand of Eq. (5.6) by the effective Fermi velocity 
vpiv) dressed with the self-energy corrections in Eq. (4.31). It can be easily seen that 
expanding the latter in powers of Cq corresponds to generating the iteration of self-energy 
corrections to the internal fermion propagators in the equation of Fig. 4. 

A solution of the form (5.7) can be still found for r3(0, 0; k, cj/;), where now each term 
in the series can be obtained from all the precedent by expanding vf{p) in powers of 
in Eq. (5.6) and using repeatedly the formula (5.8). The Fermi velocity vp needs to be 
renormalized to absorb the divergence of vpiv) iii the limit e — )■ 0, for which we define 

Vp = ZyVp^ren (5.27) 

As in subsection 4.2, just contains a simple pole 

Z^ = l + hi^ (5.28) 

with hi = —e^ /IQiTVp^vcn- After subtraction of the self-energy pole, the rest of poles in Fa 
must be reabsorbed by the multiplicative renormalization (5.9), where now Z^ is given by 
a different series 

^™ = 1 + E ^ (5.29) 

One can see that in this case Fs^rcn can be made also finite in the limit e — t- 0, with a 
set of functions di[\) that only depend on the renormalized coupling constant 

A = M"^Z„Ao = 1 (5.30) 
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The first orders in the perturbative expansion are given for instance by 

~ 1 1 2 vr2 + 1201og^(2) ^3 10^2log(2) + 6881og3(2) + 15C(3) ,4 
d,{X) = --X--log{2)X A A 

ISvr^ + 20647r2 log2(2) + 144 (7161og^(2) + 371og(2)C(3)) . 

737280 

"294^ (l69vrMog(2) + 567744 log^ (2) + 49320 log2(2)C(3) 

+57r2 (28641og^(2) + 37C(3)) + 1125C(5)) A*' + . . . (5.31) 

= 1 + 1 log(2) A» + 5^!±Iili5l!P) ,4 

^ ^ 16 24 ^ 18432 

IIOtt^ log(2) + 8592 log^(2) + 135C(3) ^5 
184320 

2937r'^ + 589447r2log2(2) + 72 (443921og'^(2) + 17791og(2)C(3)) g 
+ 53084160 ^ ^^■^^> 

d3(A) = -^log(2)A^-^^!±^^^l^^A^ 
^ 768 ^ 184320 

IOOtt^ log(2) + 11904 log3(2) + 45C(3) 

4423680 ^ ' ' 

rf5(A) = -^log(2)A6 + ... (5.35) 

We have computed the expansions of the functions dj(A) numerically up to order A^^, 
checking that the coefficients do not depend on the momentum k of the vertex. This is 
the essential requirement guaranteeing the renormalizability of the theory, by which the 
divergences in the limit e — )■ can be absorbed into the redefinition of a finite number of 
local operators. 

We can proceed to the computation of the anomalous dimension 7^ in the present 
approach, taking again into account that the scale dependence of Ta rcn stems from Zm, so 
that 

/i dX dZ,n 

^- = ^-d-^-Jx 

The dependence of the renormalized coupling A on /i can be obtained by differentiating 
(5.30), leading to 

dx X dxdz. 

Using the fact that Z„ = 1 + X{dZ^/dX), we obtain 

= -eX 

= -eA-A6i(A) (5.38) 
This result can be now introduced in Eq. (5.36), finding that 

Im = ^{-eX - X 61(A)) -^d.(A) (5.39) 

^ i=l 
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Equivalently, we can write 

f 1 + E ^1 7». = -AE ;:,;^*+i(A) - A(„(A) Y. H^dM (5.40) 



1=1 / 1=0 1=1 



Assuming that 7^ must have a finite hmit as e — 0, we get 

Im = -A-^(ii(A) (5.41) 
aX 

together with the conditions for the cancellation of all the pole contributions to jm [32] 

MX) ^di(A) = ^d,+i{X) + h{X) ^<i,(A) (5.42) 

Quite remarkably, we have verified that the conditions (5.42) are indeed satisfied by 
the perturbative series of the functions di{X), at least up to the order A^^ we have been 
able to carry out numerically the expansion. This guarantees that Eq. (5.41) can be used 
to obtain a finite result for 7^, only dependent on the value of the coupling constant A. 
Computing the perturbative expansion 



Ji(A) = ^4")a" (5.43) 

n=l 



~(n) 

we have checked that in this case again the coefficients d^ grow geometrically with the 
order n, as shown in Fig. 7. The ratio d^^~^^^ /d^^ converges to a limit value d, in such a 
way that 

CO 

- Ji(A) = ^(i"A" + regular terms (5.44) 

n=l 

We have found that the points d^^'^^^/J^"^ can be fitted quite accurately by the scaling 
behavior 

4"^ J d' d" d'" 
^ - d+ — + ^ + ^ + ... (5.45) 



We obtain in this way a finite radius of convergence for the perturbative expansion 

Ac = i w 0.544775 (5.46) 
d 

where the error is estimated to be in the last digit (in similar fashion as for the critical 
value given in (5.25). This represents a considerable gain in precision with respect to the 
calculation reported in Ref. [31], which is now made possible by our ability to carry out the 
power series expansion for di(A) up to seven orders beyond that achieved in that paper. 

The existence of a critical coupling Ac implies that the many-body theory of interacting 
Dirac fermions enters a new phase for sufficiently large strength of the Coulomb interaction. 
According to our preceding discussion, this phase is characterized by the condensation of the 
staggered charge density, which signals the dynamical breakdown of the chiral symmetry 
of the theory. We find therefore that this phenomenon takes place even after accounting 
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Figure 7. Plot of the absolute value of the coefficients d " in the expansion of di(A) as a power 
series of the renormalized coupling A. 

for the electron self-energy corrections in the many-body theory. Thus, the effect of Fermi 
velocity renormahzation tends to reduce the effective interaction strength, but it does not 
prevent the dynamical mass generation, leading instead to a larger value (5.46) of the 
critical coupling in comparison to the approach without switching on the electron self- 
energy corrections (see also Ref. [41]). 

6 Conclusions 

The above results provide strong evidence that the model of interacting Dirac fermions in 
graphene constitutes a completely renormalizable field theory, in the sense that all the cutoff 
dependences can be absorbed into the redefinition of a finite number of local operators. 
We have proven this fact for different vertex functions in the ladder approximation, as 
well as when this is supplemented by electron self-energy corrections to the electron and 
hole states in the vertices. We have seen that there is a nontrivial cancellation of poles 
in the dimensionally regularized anomalous dimensions of different fermion bilinears, a 
feature typically enforced by nonperturbative equations for the residues of the poles in 
renormalizable quantum field theories. It would be then interesting to pursue the program 
of renormahzation in other approximations with higher diagrammatic content, testing the 
cutoff independence of anomalous exponents in the graphene electron system. 

The other important conclusion we reach is that the incorporation of electron self- 
energy corrections is in general required to preserve the gauge invariance that relates the 
kinetic and the interaction terms in the effective action of the theory. When looking at 
corrections to the current density vertex in the ladder approximation, gauge invariance 
implies in particular the complete cancellation between the cutoff-dependent part of self- 
energy and vertex corrections to all orders in the perturbation expansion, as we have 
checked explicitly in our calculation with dimensional regularization. 

From a practical point of view, our computational approach has allowed us to address 
the question of the chiral symmetry breaking in the interacting theory of Dirac fermions in 
graphene. Making use of the renormalizability of the theory, we have shown that quantities 
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like the anomalous dimension of the staggered charge density operator can be determined in 
terms of the renormalized coupling alone, allowing to characterize the dynamical symmetry 
breaking at the critical coupling given by the finite radius of convergence of the perturbative 
expansion. 

We note anyhow that the critical value found for the coupling A cannot be used directly 
to predict the onset of dynamical symmetry breaking in real graphene, as the e-e interaction 
is affected in general by screening processes that may reduce significantly its strength. This 
means that the critical coupling Ac we have computed should be referred actually to the 
effective interaction strength after incorporating screening corrections. These depend in 
general on intrinsic factors, like the number N of different fermion flavors that enter in the 
polarization of the system. Under static RPA screening, for instance, the coupling A of 
the effective interaction is related to the bare fine structure constant a in graphene by the 
expression 

A = ^ — 6.1) 

N being the number of four-component Dirac fermions. In the physical case N = 2, the 
nominal coupling of graphene in vacuum a ~ 2.2 leads to the estimate A ~ 0.49, which 
is above the critical coupling (5.25) obtained in the pure ladder approximation, but below 
the value (5.46) found after incorporating the electron self-energy corrections. The latter 
critical value of A translates therefore into a more stringent bound on the nominal value 
of jvp for the development of chiral symmetry breaking. It is worth noticing however 
that the more sensible screening approaches, considering the dynamical polarization of 
electron-hole pairs, still predict the dynamical mass generation for the coupling a ~ 2.2 
corresponding to graphene isolated in vacuum[29, 31]. 

The reliability of our renormalization procedure is reinforced by the fact that the 
critical coupling we obtain from the sum of ladder diagrams for the vertex matches with 
great accuracy the value found within a quite different approach to chiral symmetry break- 
ing in graphene, based on the self-consistent resolution of the gap equation for the Dirac 
fermions [27]. This provides good motivation to check whether a similar agreement can be 
reached after incorporating the electron self-energy corrections in the self-consistent gap 
equation, as well as to extend by any other feasible means the approach devised in the 
present paper. 
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